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ABSTRACT 


A  linearly  polarized  antenna  cannot  radiate  power  uniformly  in 
all  directions.  However,  by  controlling  the  aperture  excitation,  as 
is  done  in  an  array,  it  is  possible  to  reduce  the  maximum  gain  and 
approach  conditions  for  isotropic  radiation.  A  numerical  method  is 
presented  which  generates  a  family  of  designs  which  depend  on  a  param¬ 
eter  a.  As  a  approaches  zero,  the  radiation  pattern  tends  to  become 
more  isotropic.  However,  the  efficiency  is  reduced  and  the  sensi¬ 
tivity  of  the  pattern  to  errors  in  the  aperture  function  is  increased. 
In  some  cases  this  sensitivity  is  so  high  as  to  make  the  result  worth¬ 
less.  The  design,  therefore,  must  be  a  compromise  between  closeness 
to  conditions  of  isotropic  radiation  on  one  hand,  and  losses  and  sen¬ 
sitivity  to  errors  on  the  other. 

The  sensitivity  to  errors,  expressed  by  the  pattern  deteriora¬ 
tion  for  a  given  level  of  error,  has  been  evaluated  by  a  numerical 
experiment  (Monte  Carlo  method).  A  general  relation  has  also  been 
established  between  the  sensitivity  and  the  losses  in  the  antenna. 

It  agrees  with  the  numerical  experiment  and  can  be  used  as  a  guide 
to  choose  the  regularization  parameter  a. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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1-  INTRODUCTION 

A  linearly  polarized  antenna  cannot  radiate  power  uniformly  in  all 
directions.  This  condition  of  isotropic  radiation,  or  unity  gain,  is 
also  impossible  if  one  specifies  the  axial  ratio  of  the  polarization 
ellipse. ^  However,  accepting  the  existence  of  nulls  in  the  radiation 
pattern,  one  may  ask  if  it  is  possible  to  reduce  the  total  solid  angle 
in  which  the  gain  falls  below  a  certain  level.  This  is  equivalent  to 
reducing  the  maximum  gain. 

It  will  be  shown  that  it  is,  indeed,  possible  to  make  the  maximum 
gain  as  close  to  unity  as  one  wishes.  However,  the  practical  difficul¬ 
ties  in  doing  this  are  the  same  as  those  encountered  in  trying  to  increase 
the  gain  of  an  antenna  indefinitely .  To  reach  a  certain  level  of  per¬ 
formance  (gain  or  lack  of  gain),  one  must  have  an  antenna  of  sufficient 
aperture.  If,  for  a  given  aperture,  we  try  to  improve  the  performance 
above  a  certain  level,  we  have  to  increase  drastically  the  antenna  cur¬ 
rents  (or  the  near  field).  As  a  consequence,  both  the  losses  and  the 
sensitivity  to  construction  errors  are  increased.  This  is  character¬ 
istic  of  the  supergain  problem.  It  occurs  more  generally  anytime  one 
tries  to  shape  a  radiation  pattern  too  accurately.  In  particular,  it 
occurs  in  the  ’’minimum  gain"  problem  considered  here. 

The  need  for  an  isotropic  radiator  arises  when  one  disposes  of  a 
certain  amount  of  power  and  one  wishes  to  produce  a  radiation  intensity 
higher  than  a  given  level  in  a  solid  angle  as  large  as  possible.  This 
may  be  for  the  purpose  of  reaching  a  target  whose  direction  is  not  known 
a  priori.  Conversely,  in  the  receiving  case, one  may  want  to  receive  a 
sufficient  signal  from  a  source  of  given  power  and  unknown  direction. 
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Although  the  method  used  could  be  generalized,  we  shall  consider 

in  this  report  a  colinear  array  of  dipoles  as  shown  in  Figure  la.  From 

the  far-field  power  pattern  of  this  antenna,  which  is  shown  in  Figure  lb, 

and  some  gain  level  G  ,  one  can  read  the  total  solid  angle  ft  in  which 

m 

the  gain  falls  below  G  .  The  plot  of  G(ft)  versus  ft  is  shown  in  Figure 

m 

lc.  In  this  plot,  (4tt— ft)  is  the  total  solid  angle  in  which  the  gain 

exceeds  G  . 

m 

Given  the  curve  G(ft),  the  total  available  power  P  watts,  and  the 
desired  radiation  intensity  I  watts/st eradian,  the  useful  angle  (47T-ft) 
is  defined  by 


G(4ir-n)  =  I 

or 

47tI 

G(4ir-ft)  -  =  G  . 

r  m 

If  the  ratio  I/P  is  small,  ft  will  be  small  for  any  design.  However, 
when  I/P  is  large,  then  special  designs  become  necessary  in  order  to 
produce  an  intensity  I  in  as  large  a  solid  angle  as  possible. 

In  Figure  lc,  the  shaded  portion  under  the  gain  curve  has  an  area 
equal  to  4tt.  The  total  area  of  the  rectangles  ABCD  and  DEFH  exceeds 
that  of  the  shaded  portion.  Accordingly,  the  following  inequality  may 
be  written, 


( 4 7T — ft )  Gw  +  olG  >  4tt 
M  m  — 

or 


ft  1  4*  (Gn-1)/  (Gm-g)» 
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(C) 


Figure  1.  Schematic  representation  of  array  and  its  far-field  gain 
and  power  patterns. 

(a)  Array  of  colinear  dipoles. 

(b)  Far-field  power  pattern. 

(c)  Curve  of  gain  versus  solid  angle. 
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where  is  the  maximum  gain.  This  last  inequality  implies  that  for 
ft  0,  the  maximum  gain  must  approach  unity. 

To  obtain  an  indication  of  what  can  be  achieved  practically  with  a 
moderately  short  antenna  (1  wavelength) ,  this  report  considers  the  syn¬ 
thesis  of  linear  antenna  arrays  made  of  dipoles  aligned  in  some  direc¬ 
tion,  Ox,  as  shown  in  Figure  la.  The  resulting  pattern  has  nulls  in  the 

directions  +x.  The  problem  is  to  reduce  the  cones  where  the  gain  is 

below  a  certain  level.  A  general  method  based  on  the  concept  of  regular- 
(2  3) 

ization  *  was  developed  in  the  Antenna  Laboratory  under  NSF  Grant 
GK  2120.  It  will  be  applied  to  the  problem. 
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2.  STATEMENT  OF  THE  PROBLEM 

The  general  case  of  a  pattern  function  g^(0  defined  over  the  vis¬ 
ible  region[-k  _<  £  <_  k]  will  be  first  -considered,  and  the  numerical 
results  will  be  specialized  to  the  case 

gQ(0  =  1  [-k<Elk]  (1) 

where  the  gain  is  equal  to  one  over  the  whole  visible  region.  A  source 
distribution  f(x)  representing  an  array  of  N  elements  over  the  aperture 
[-a  _<^  x  <_  a]  is  of  the  form 

N 

f(x)  =  £  f  6 (x-x  ).  (2) 

,  n  n 
n^l 

The  problem  is  to  determine  the  {f^}  such  that  the  pattern  function 
of  the  array  comes  close  to  gg(£). 

The  far-field  pattern  of  a  linearly  polarized  source  f(x)  is  given 
by 

a  i£x 

g(5)  =  Tf  =  e(£)  /  f (x)e  dx  [-k  _<  £  _<  k]  (3) 

-a 

where  £  =  k  sin  9.  The  first  factor  e(£)  is  the  element  pattern. 

For  a  source  distribution  described  in  Equation  (2) ,  the  far-field 
pattern  is  given  by, 

N 

g(5)  -  Tf  =  Z  f  g  (o  (4) 

i  n  n 


where 

g  (O  =  e(C)ei^Xn  [-k  <  £  <  k]. 
n  —  — 
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The  function  g(£)>  therefore,  belongs  to  a  finite  dimensional  vector 
space  0  spanned  by  the  N  elements  g^CO  (1  <  n  <  H).  In  general,  &q(0 
does  not  belong  to  this  subspace  and  the  problem  is  to  determine  the 
complex  coefficients  { fn }  such  that  the  pattern  given  by  (4)  comes  close 
to  gQ(C) • 

In  one  method  proposed  in  the  next  Chapter,  the  coefficients  {fn} 

are  determined  such  that  the  corresponding  pattern  g(C)  is  the  orthogonal 

N 

projection  of  &q(0  on  the  finite  dimensional  subspace  0  .  This  gives 
the  least-squares  approximation  to  • 

This  straightforward  projection,  however,  may  result  ip  a  solution 
with  a  high  Q  factor,  that  is,  excitation  currents  { f that  are  very 
large  for  a  given  power  radiated.  These  currents  vary  widely  from  ele¬ 
ment  to  element,  and  the  resulting  far-field  pattern  becomes  very  sen¬ 
sitive  to  errors  in  construction. 

(2) 

To  avoid  this  difficulty  in  problems  of  the  same  type,  Tihonov 

has  proposed  a  method  of  regularization  which  consists  of  imposing  some 

restrictions  on  the  acceptable  solutions,  that  is, on  the  aperture  functions  f. 

The  regularization  algorithm  depends  on  a  parameter  a  such  that  for 

N 

a  =  0,  we  have  the  direct  projection  on  the  subspace  t  •  As  a  increases, 
a  set  of  solutions  is  obtained  with  decreasing  values  of  the  Q  factor. 

This  set  of  solutions  invariably  produces  far-field  patterns  that  are 
further  away  from  the  desired  pattern  but  are  also  less  subject  to 
ohmic  losses  and  less  sensitive  to  errors  in  construction. 

To  give  a  precise  meaning  to  the  concept  of  projection  and  of  Q 
factor,  one  has  to  introduce  appropriate  norms  in  the  space  3~  of  source 
functions  f  =  {f^}  and  the  space  (Y  of  far-field  patterns  g(£).  In  both 
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spaces,  it  is  natural  to  choose  the  L ^  norm  which  can  be  interpreted 
in  the  space  as  the  power  radiated  by  the  antenna.  Thus, 


N 

E 

n=l 


f  f 
n  n 


(5) 


and 


I  I  si  I2  =  7  /  g*(C)  g(0  dg.  (6) 

k  -k 

The  Q  factor  for  a  given  design  is  represented  by 

Q  =  1 1 f 1 12  /  ||g||2.  (7) 

Because  of  the  element  pattern,  this  factor  is  not  always  larger  than 
one. 
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3.  METHOD  OF  SOLUTION 


Tihonov's  procedure  consists  of  minimizing  the  functional 
=  |  |g-gQ|  I2  +  a|  |f |  |2 


(8) 


where  a  is  a  positive  constant  referred  to  as  the  regularization  param¬ 
eter. 


By  setting  the  first  variation  of  J^(f;g)  with  respect  to  f  ■  {f^} 

equal  to  zero,  the  following  equation  for  f  =  {f  },  the  minimizing  set 

n 

of  excitations,  is  obtained, 


N  ,  k  * 

If  f/ggdC+af 
,0k.  &n°m  a 

n=l  n  -k  m 


k  {  80  ’  8m  d5 
-k 


(1  <  n ,  m  <  N) 


where,  since  is,  symmetrical ,  the  excitations  {f^}  are  assumed  real. 
In  matrix  notation, 

(G+cI)fo  -  G0g0  (9) 

where  G  is  an  (N  x  N)  square  matrix  with  terms  G  given  by 

1  k  *  i£(x  "X  ) 

G  -  -r  I  e(Oe(C)e  n  d5.,  (10) 

mn  k  . 

-k 

and  I  is  the  identity  matrix  of  dimension  N. 

The  column  matrix  g^  =  (gg  }  (1  £  J  i  represents  the  sam- 

j 

pled  values  of  the  function  g^  at  the  M  equally  spaced  points  over  the 
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visible  region  [-k  <_  £  _<  k],  and  Gq  is  a  (N  x  M)  quadrature  matrix  derived 
to  perform  numerically  the  integral 


1 

k 


k 


f  g0gm*d? 


for  any  given  function  g^(C) • 

Equation  (10)  shows  that  the  G  are  the  Fourier  Transforms  H(x)  of 

mn 

2 

h(£)  =  |e(Ol  evaluated  for  x  -  x  . 

n  m 


G 

mn 


H(x  -x  ) 
n  m 


When  the  elements  are  dipoles  in  the  x  direction, 

2  1/2 

e(0  =  cos  0  =  (1 - j  ) 

k 

and 

J-,  (kx) 

H(x)  =  2/(2*)  J 

(kx) 


(ID 


4  r  sin(kx) 
(kx)2  (kx) 


cos (kx) ] . 


(12) 


To  derive  the  elements  Gn  of  the  Quadrature  matrix  Gn,  the  func- 

ij 

tion  is  sampled  at  equally  spaced  points  over  the  visible  region 
[k  _<  £  _<  k]  and  approximated  by  straight  lines  over  each  interval.  For 
the  range , 


cj-i  - 5  -  5j 


g0(s)  •  £  1(c-c)l 
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where  h  is  the  length  of  each  interval.  For  the  range, 


~  ^  1  Cj+1 


s0(O  h  [sj+i(c"cj)_8j(c‘S+i)3 


and  the  terms  are  given  by 

ij 


1  1 


J  2/2 

°0.  'th</  <«-Vl,<1'f2)  6  « 

1J  S-l  , 

j+1  ^2  \/2  ^xi 

-/  (c-c,+1)d  -  %  )  e  d?}. 

k 

Since  these  integrals  cannot  be  evaluated  in  closed  form,  they  were  eval¬ 
uated  by  Simpson's  three-point  quadrature.  From  (9), 


-1 


f,  -  (G+oI)  G0g0 


(13) 


Once  a  choice  for  a  is  made,  Equation  (13)  gives  a  unique  solution 
f^  for  a  given  g^.  This  solution  f^  has  a  far-field  pattern  g^,  given 
by 

g  =  Tf  (14) 

a  a 

and  a  corresponding  Q  factor.  The  function  f^  has  the  property  that, 
for  this  value  of  Q, its  far-field  pattern  yields  the  least-squares 
fit  to  gQ. 


4.  NUMERICAL  RESULTS 
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In  order  to  study  the  relationship  of  ot  to  the  sensitivity  of  the 
design  to  errors,  Equation  (13)  was  solved  for  a  number  of  values  of  a 
ranging  from  0  to  10  \ 

The  magnitudes  of  the  element  excitations  for  arrays  of  five  and 
seven  elements  are  reported  in  Tables  I  and  II,  respectively ,  for  various 
values  of  the  regularization  parameter.  The  element  pattern  e(£)  is 
given  in  Equation  (11). 

TABLE  I.  'MAGNITUDE  OF  ELEMENT  EXCITATIONS  FOR  THE  FIVE-ELEMENT  ARRAY. 


Element 

Number 

a  =  0 

a  =  10  3 

in"2 

a  =  10  a 

ii 

i— * 

o 

i 

H* 

i 

.64 

.35 

-.15 

-.28 

2 

-2.25 

-1.47 

-.11 

.335 

3 

4.28 

3.26 

1.47 

.86 

4 

-2.25 

-1.47 

-.11 

.336 

5 

,64 

.35 

-.15 

-.28 

The  far 

uated  using 

-field  pattern 

Equation  (4) . 

g^  corresponding  to 

Since  it  is  assumed 

each  solution  f 

ot 

that  only  power 

is  eval- 

2 

P  =  |  Ig^l  |  is  available  to  the  array  for  radiation,  the  excitations 

{f  }  are  normalized  such  that 

ot 

n 

Let  f  =  {f  }  and  g  be  the  normalized  sources  and  far-field  patterns, 
a  a  ot 

n  2 

respectively.  The  relative  error  due  to  regularization  is  expressed 

_  2 

by  the  norm  of  the  error  | | g^— gg | |  divided  by  that  of  the  desired 
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TABLE  II.  MAGNITUDE  OF  ELEMENT  EXCITATION  FOR  THE  SEVEN- ELEMENT  ARRAY 


Element 

Number 

a  =  0 

a  =  10  6 

a  =  10  5 

P 

1 

I-1 

o 

1 

-p> 

i 

-6.71 

-1.58 

.57 

.87 

2 

35.07 

10.19 

-.29 

-1.80 

3 

-81.46 

-26.53 

-3.36 

.07 

4 

107.20 

36.88 

7.22 

2.79 

5 

-81.46 

-26.53 

-3.36 

.07 

6 

35.07 

10.19 

-.29 

-1.80 

7 

-6.71 

-1.58 

.57 

.87 

Element 

Number 

a  «  10  2 

a  =  10_1 

1 

.05 

-.23 

2 

-.55 

-.01 

3 

.43 

.41 

4 

1.08 

.62 

5 

.43 

.41 

6 

.55 

-.01 

7 

.05 

-.23 
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Next,  in  order  to  study  the  sensitivity  of  each  design  to  errors  in 
construction,  a  random  relative  error  v_  is  introduced  to  the  source  func¬ 
tion,  and  these  excitations  with  noise  are  denoted  by  f  and  given  by 

n 

f  “  f  (1+v  )  (16) 

av  a  n 

n  n 

where  _v  =  {v  }  are  random  complex  numbers  with  a  Gaussian  distribution 


E(v  )  -  0 
n 

*  2 
E(v  v  )  -  v  6  . 

n  m  nm 

The  far-field  pattern  corresponding  to  sources  {f  }  is  denoted  by  g 

av  av 

n  2  2 

The  excitations  {f  }  are  normalized  such  that  Mg  II  =  !  1  gA  I  . 

av  1  av 11  1 1  O'  1 

n_ 

Let  f  ■  {f  }  and  g  be  the  normalized  sources  and  far-field  pattern, 

av  av  av 

n  2 

respectively.  The  relative  error  e  ,  due  to  the  introcjuction  of  random 


noise,  is  expressed  by  the  norm  of  the  error  | !g  -g^| 


divided  by  that 


of  the  desired  pattern. 


\2  ■  H8ou-iJI2  / 


(17) 


and  the  corresponding  relative  errors  due  to  regularization  and  the  intro- 

^  2 

duction  of  random  noise  are  expressed  by  the  norm  of  the  error  | |g^  -g^ j  | 
divided  by  |  | gQ | | 2: 


2 


e 


t 


/  I  I  Bq 


2 


For  each  design,  the  Q-factor  is  calculated  with  the  presence  of  random 

errors  and  denoted  by  Q  ,  where 

v 


Q  = 


(18) 


14 

In  each  instance,  25  sets  of  random  errors  were  introduced  and  the  aver- 


2  2 

age  of  the  quantities  and  (1/Q  )  were  calculated  over  the  25  trials, 

where 


2  2 

e  =  <e  > 
v  v  av 


2 


e 

t 


2 

<e  > 

t  av 


(19) 


1/Q  =  <1/ Q  >  . 

v  av 

The  quantities  e^,  e  ,  and  are  illustrated  schematically  in  Figure  2. 

In  this  diagram,  vectors  are  understood  to  represent  functions  and  straight 

lines,  distances  in  the  sense  of  the  norm.  The  curved  lines  represent 

operators  mapping  the  space  of  source  functions  into  the  space  of  far- 

field  pattern  functions.  The  three  vectors  gn,  g  ,  and  g  are  drawn 

0  a  av 

with  their  tips  on  the  same  circle  to  indicate  that  their  norms  are 
equal. 

The  far-field  power  patterns  ||g  ||  and  |  g^  |  of  the  arrays  for 
various  values  of  the  regularization  parameter  a  are  shown  in  Figures  3  to  11. 
Also  in  the  same  figures,  the  power  pattern  intensity  of  an  elementary 
dipole  is  shown  for  comparison.  The  horizontal  line  drawn  for  a  radia¬ 
tion  intensity  «  1  is  that  of  an  ideal  isotropic  radiator.  It  is 
pointed  out  again  that  all  pattern  functions  are  normalized  such  that 
their  norms  are  equal.  This  makes  their  radiated  power,  represented 

by  the  area  under  the  curve,  equal  in  all  cases.  In  Figures  12  and  13, 

2  2  ~~ 2 

curves  for  the  quantities  e  ,  e  ,  e  and  Q  are  shown  as  a  function  of  the 

a  v  t 

regularization  parameter  a  for  the  arrays  of  five  elements  and  seven  ele¬ 
ments  ,  respectively . 


Space  of  Pattern  Functions  g 
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Figure  2.  Schematic  representation  of  the  quantities 
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A1ISN31NI  U3M0d  Q1 3I3-UV3 


SOLID  ANGLE  CC  £  sin  9 
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A1ISN31NI  d3M0d  Q1 313-dVd 
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A1ISN31NI  d3M0d  Q1 313-UVd 
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SOLID  ANGLE  CC  £  =  sin  6 
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REGULARIZATION  a 


Figure  12. 


Error  parameters  e^,  ev,  and  e^,  and  Q  vs.  regularization  parameter 
a  for  the  five-element  array. 
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Figure  13. 


Error  parameters  ev,  and  and  Q  vs.  regularization  parameter 
a  for  the  seven-element  array. 
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2 

From  the  far-field  power  pattern  intensity  |g^|  in  Figures  3  to  11, 

curves  for  the  gain  vs.  solid  angle,  as  described  in  the  introduction, 

were  derived.  These  curves  are  shown  in  Figures  14  and  15  for  the  five1- 

-2 

and  seven-element  array,  respectively,  for  the  cases  a  =  0  and  a  -  10 
The  gain  curve  for  an  elementary  dipole  is  also  shown  for  comparison. 

For  example,  in  Figure  14,  for  a  directive  gain  level  of  0.8,  the  total 
solid  angle  in  which  the  gain  of  each  of  the  cases  considered  exceeds 
the  0.8  gain  level  is  as  follows: 


Ideal  isotropic  radiator  : 

4n 

Five  element  array 

(a  =  0)  : 

3.82tt 

Five  element  array 

(a  =  10  : 

3. 36tt 

Elementary  dipole 

2. 72it. 

A  plot  of  1/Q  vs.  1/Q  is  shown  in  Figure  16. 


Element  Pattern 
(cos  6) 
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Figure  14.  Gain  vs.  solid  angle  for  the  five-element  array. 


Element  Pattern  (cos 9) 
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UID9  aAipajjQ 


Figure  15.  Gain  vs.  solid  angle  for  the  seven-element  array. 
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Figure  16.  Plot  of  1/Qy  vs.  1/Q  for  both  five-  and  seven-element  arrays 
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5.  COMMENTS  ON  THE  NUMERICAL  RESULTS  AND  THE  CHOICE  OF  a 

The  first  observation  is  that  using  an  array  does, indeed,  improve 
the  low  gain  performance.  The  gain  curves  in  Figures  14  and  15  for 
a  =  0  come  closest  to  that  of  an  ideal  isotropic  radiator  among  all 
designs  considered.  Comparing  the  gain  curves  of  the  five-element  to 
the  seven-element  array,  only  a  slight  improvement  in  low-gain  perform¬ 
ance  is  observed  in  the  latter.  The  aperture  size  in  both  cases,  how¬ 
ever,  is  the  same. 

The  performance  of  those  designs  corresponding  to  relatively  small 

values  of  a  looks,  indeed,  promising.  However,  in  the  presence  of  random 

errors  in  the  currents,  the  performance  of  these  designs  deteriorates 

drastically,  as  shown  in  Figures  3  to  4  and  6  to  9.  In  such  instances, 

a  single  element  has  a  better  performance  than  the  array! 

In  deciding  which  value  of  a  to  use,  it  becomes  imperative  to  study 

2  2 

the  behavior  of  the  quantities  e^,  and  Q.  For  a  constant  power  radi¬ 

ated,  the  Q-f actor  is  directly  proportional  to  ohmic  losses  while  the 
2 

error  is  an  indication  of  the  sensitivity  of  the  design  to  error. 

2 

In  Figures  12  and  13,  as  a  approaches  zero,  decreases,  and  both  Q  and 
2 

increase.  For  the  larger  values  of  a,  these  trends  are  reversed. 

2  2 

The  design  then  must  be  a  compromise  between  e  ,  and  Q.  A  good  com- 

— 2 

promise  is  that  value  of  a  which  minimizes  the  total  error  For  the 

_3 

five-element  array,  this  corresponds  to  a  =  3.54  x  10  ,  Q  =  2.25,  and 

2  -3 

=  .008  and,  for  the  seven-element  array,  to  a  =  1.6  x  10  ,  Q  =  4.0, 

2 

and  =  .0045.  For  these  particular  values  of  a,  the  seven-element 
array  has  larger  ohmic  losses,  but  is  less  sensitive  to  errors  than  the 
five-element  array . 
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It  is  observed  in  Figure  16  that  the  introduction  of  random  errors 

in  the  element  excitations  produces  a  decrease  of  the  Q-factor,  that 

is,  an  increase  of  the  power  radiated  for  a  given  norm  of  the  sources. 

For  values  of  Q  less  than  10,  the  reduction  is  very  small.  For  values 

of  Q  larger  than  10,  the  reduction  becomes  appreciable , and  for  values 

4 

of  Q  larger  than  10  Q^, becomes  cons tant,  equal  to  270. 
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6.  STATISTICAL  ESTIMATION  OF  1/Q  AND  e2 

V  V 

2 

The  numerical  determination  of  1/Q  and  e  ,  as  outlined  in  the  pre- 

v  v 

vious  section,  can  be  time  consuming  since  the  averages  have  to  be  taken 

over  a  large  number  of  trials.  It  would  be  advantageous,  therefore,  if 

2 

the  quantities  1/Q  and  could  be  estimated  analytically.  From  the 

expression  for  f  in  Equation  (16)  and  the  definition  of  the  L.  norm 

av  2 

in  the  space  of  aperture  source  functions  given  in  Equation  (5) ,  the 

norm  of  f  is  evaluated: 

av 


N 

E 

n=l 


f  f 

a  .  a 


n  n 


(l+v  )(l+v  *), 
n  n 


and  the  expected  value  of  this  quantity  is  found  to  be 


E{ 


If  || 

av 


|f  | I2  (l+v2) 
'  a  1  1 


where, for  v  <<  1 


E{  |  |  f 


av 


(20) 


The  far-field  pattern  g  ,  corresponding  to  an  aperture  source  distribu¬ 
tion  f^,is  given  by  Equation  (4): 


av 


=  "  E  f 
n=l 


n 


From  Equations  (6)  and  (10),  the  norm  of  g  is  evaluated, 


N  N 
E  E 

n=l  ra=l 


-  * 
G  f 
nm  a 


* 


)  a+vm> , 
m 


and  the  expectation  is  found  to  be 


wllsjl2)  ■  IliJI2  +  v2g00i |fj I2. 


(21) 


The  expected  value  of  1/Q  is  defined  to  be 


E{l/Q  }  =  E{ | | g  | | 2  /  | | f  ||2}. 
v  1  1  av  1  1  1  av 1  1 


However,  Equation  (20)  shows  that  when  the  mean  square  level  of  error 

2  “2 

v  is  small,  the  expected  value  | | f  | |  may  be  approximated  by  | |f  | | 
and  the  following  relation  for  1/Q^  is  obtained  from  Equation  (21) : 

1/QV  =  E{1/Qv>  - l/Q  +  Gqqv 2 .  (22) 

2 

Before  computing  the  relative  error  e  due  to  random  errors,  g 

v  °ctv 

should  be  normalized  so  that 


El  |'|g 


av  1 


2 


From  Equation  (21) ,  the  normalization  factor  F  is  given  by 


F  =  [1+v2G0()Q]  1/2  =  [Q-E  {  ~  }  ]  1/2.  (23) 

v_ 

Let  g  denote  the  normalized  far-field  pattern  of  the  aperture  sources 
in  the  presence  of  random  error. 


6cv  =  F  6av- 

2 

The  relative  error  is  given  by 
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Figure  17.  Schematic  representation  of  the  relationship  between  the  pattern 
functions  and  the  error  parameters. 
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Figure  18.  Plot  of  the  theoretical  and  numerical  values  of  e  and  1/Q 
vs.  the  Q-factor. 
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where 


N 

g  -  g  “  E  f  g  [ F ( 1+v  )-l]. 
av  a  .  a  &n  n 

n=l  n 

Taking  the  norm  and  then  the  expectation,  the  following  relationship 
2 

is  obtained  for  e  : 

v 

ej  =  2 (1-F) .  (24) 

2 

This  equation  relates  simply  to  the  drop  in  the  Q-factor  when  random 
errors  are  introduced  in  the  sources. 

Equation  (24)  can  be  rederived  from  the  schematic  representation 
in  Figure  17.  Equation  (21)  shows  that  the  quantities  ||g^||,  v |  | f  |  |  /Gqq , 
and  [E{||g  ||  }]  are  related  by  the  Pythagorean  theorem.  Therefore, 
the  tangent  of  the  angle  0  shown  in  Figure  17  is  given  by 

tan  6  =  vv'qGqq  (25) 

and 

e  ■  2  sin  —  ,  (26) 

v  2 

which  can  be  rewritten  in  a  form  similar  to  that  of  Equation  (24) : 

=  2(l-cos  0). 

v 

2 

Plots  of  1/Q^  and  e  as  a  function  of  Q  calculated  from  Equations 
(22)  and  (24)  are  shown  in  Figure  18.  The  values  of  these  quantities 
evaluated  numerically  are  also  shown  for  comparison.  The  agreement  is 


very  good. 
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7.  CONCLUSION 

A  method  has  been  outlined  to  design  a  linear  array  whose  gain  is 
close  to  unity  and,  therefore,  whose  pattern  is  nearly  isotropic.  Tihonov’s 
regularization  algorithm  is  used  to  produce  a  family  of  designs  which 
depend  on  a  parameter  a. 

2 

It  is  found  that  as  a  approaches  0,  the  error  from  the  actual  to 

2 

the  desired  pattern  decreases.  But  the  norm  |f||  of  the  aperture 

function,  which  for  a  unit  power  radiated  is  also  the  supergain  ratio 

Q  introduced  by  Taylor,  increases.  This  is  accompanied  by  an  increase 

2 

in  the  sensitivity  of  the  design  to  errors  represented  by  .  An  accept¬ 
able  design  has  to  be  a  compromise  between  these  three  quantities. 

Although  the  sensitivity  to  error  varies  in  the  same  direction  as 
Q,  it  is  not  directly  proportional  to  it.  A  general  relation  has  been 
established  between  these  two  quantities.  This  relation  has  been  veri¬ 
fied  numerically  and  it  provides  a  useful  guide  for  choosing  the  amount 
of  regularization . 

The  regularization  method  proves  to  be  effective  in  designing  nearly 
isotropic  arrays  of  reasonable  sizes. 
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